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Abstract 



We present a new derivation of the cosmic microwave background anisotropy 
spectrum from the Ostriker-Vishniac effect for an open, flat, or closed Uni- 
verse, and calculate the anisotropy expected in cold dark-matter (CDM) mod- 
els. We provide simple semi-analytic fitting formulas for the Vishniac power 
spectrum that can be used to evaluate the expected anisotropy in CDM mod- 
els for any arbitrary ionization history. In a flat Universe, CDM models 
normalized to cluster abundances produce rms temperature anisotropies of 
0.8-2.4 fiK on arcminute angular scales for a constant ionization fraction of 
unity, whereas an ionization fraction of 0.2 yields rms anisotropies of 0.3-0.8 
/xK. In an open and/or high-baryon-density Universe, the level of anisotropy 
is somewhat higher. The signal in some of these models may be detectable 
with planned interferometry experiments. The damping of the acoustic peaks 
in the primary-anisotropy spectrum at degree angular scales depends primar- 
ily on the optical depth and only secondarily on the epoch of reionization. 
On the other hand, the amplitude of Ostriker-Vishniac anisotropies depends 
sensitively on the epoch of reionization. Therefore, when combined with the 
estimate of the reionization optical depth provided by maps of degree-scale 



*This work is dedicated to the memory of David Schramm. 
^Email address: jaffe@physics.berkeley.edu 
^Email address: kamion@phys.columbia.edu 



1 



anisotropies, the Ostriker-Vishniac effect can provide a unique probe of the 
epoch of reionization. 

98.70.Vc,98.80.Cq 
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I. INTRODUCTION 



Standard cold dark matter (CDM) and its variants are the current leading models for 
the origin of large-scale structure. The canonical CDM model hypothesizes a primordial 
scale-free spectrum of primordial adiabatic perturbations and a critical density of cold dark 
matter. Its variants include low-density (e.g., Qq — 0.3 — 0.4) models, either in an open 
Universe or in a flat Universe with a cosmo logical constant, tilted models in which the 
power-law index n of the primordial power spectrum differs slightly from scale-free (n = 1), 
low-Hubble-constant models, or mixed-dark-matter models in which the Universe has a 
critical density but roughly 30% is in the form of hot dark matter. 

Although most of the matter in these models does not undergo gravitational collapse until 
relatively late in the history of the Universe, some small fraction of the mass is expected to 
collapse at early times. Ultraviolet photons released by this early generation of star and/or 
galaxy formation will partially reionize the Universe, and these ionized electrons will re- 
scatter at least some cosmic microwave background (CMB) photons after recombination at 
a redshift of z ~ 1100. Theoretical uncertainties in the process of star formation and the 
resulting ionization make precise predictions of the ionization history difficult. Constraints 
to the shape of the CMB blackbody spectrum and detection of CMB anisotropy at degree 
angular scales suggest that if reionization occurred, the fraction of CMB photons that re- 
scattered is small. Still, estimates show that even if small, at least some reionization is 
expected in CDM models |T]-§[: for example, the most careful recent calculations suggest a 
fraction r r ~ 0.1 of CMB photons were re-scattered 

Scattering of CMB photons from ionized clouds will lead to anisotropics at arcminute 
separations below the Silk-damping scale (the Ostriker-Vishniac effect) |§-1C]. These 



anisotropies arise at a higher order in perturbation theory and are therefore not included 
in the usual Boltzmann calculations of CMB anisotropy spectra. The level of anisotropy is 
expected to be small and it has so far eluded detection. However, these anisotropies may 
be observable with forthcoming CMB interferometry experiments [|TTJ that probe the CMB 
power spectrum at arcminute scales. 

In this paper, we present a new derivation of the anisotropy spectrum from the Ostriker- 
Vishniac effect. We then calculate the anisotropy expected in CDM models. We provide 
numerical results for the anisotropy predicted by a variety of CDM models with several 
plausible ionization histories. We also provide simple semi-analytic fitting formulas for the 
Ostriker-Vishniac power spectrum that can be used to evaluate the expected anisotropy in 
CDM models for any arbitrary ionization history. We estimate the anisotropy amplitude 
that may be detectable with forthcoming interferometry experiments. We find that the 
signal should be detectable in a number of CDM models with plausible ionization histories. 

To a good approximation, the damping of the acoustic peaks in the primary (linear- 
theory) CMB spectrum from reionization is the same for any ionization history that produces 
the same optical depth to the surface of last scatter (see, e.g., Ref. [l|]). Therefore, although 
MAP and Planck \Y2\ should provide a good estimate of this optical depth [13], they will 
not strongly constrain the ionization history, that is, whether reionization occurred earlier 
or later. On the other hand, our results indicate that the amplitude of the Ostriker-Vishniac 
anisotropy depends primarily on the ionization fraction and only secondarily on the optical 
depth. Therefore, when combined with the determination of the optical depth from MAP 
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and Planck, the Ostriker-Vishniac effect can determine the epoch of reionization. This will be 
important for precise determination of the baryon density, Hubble constant, and cosmo logical 
constant from the higher acoustic peaks in the CMB spectrum. It will also provide a unique 
window to the epoch at which structures first undergo gravitational collapse in the Universe. 



II. PRELIMINARIES 

The scale factor of the Universe, a(t), satisfies the Friedmann equations, 

- = h e(z) = # \M)(i + zf + n A + (i - n - O0(i + z)\ (i) 



- = H 2 [n A -n (i + z) 3 /2], (2) 

d 

where Hq = 100 h km sec -1 Mpc" 1 is the Hubble constant, Qq is the current nonrelativistic- 
matter density in units of the critical density, Q\ is the current contribution of the cosmo- 
logical constant to the critical density, and the overdot denotes a derivative with respect to 
time. If Qq + — 1) the Universe is flat; if it is greater (less) than unity the Universe is 
closed (open). An Einstein-de Sitter Universe has fi = 1 and il A = 0, so E(z) = (1 + zf"! 2 . 

We choose the scale factor such that a^H^ = 2. If we are located at the origin, w = 0, 
then an object at redshift z is at a comoving distance, 

1 f 2 dz' 

We define a conformal time by dr] = dt/a, so the comoving distance to the horizon is the 
conformal time today, r/ = io(oo). In an Einstein-de Sitter Universe, a/a = rf and r] = 1 
with our conventions. 

If the density contrast at comoving position w at time t is 5(w, t), with Fourier transform 

S(k) = J d 3 w expHk • w) 6(w) 5(w) = J exp(zk- w)<5(k), (4) 

then the power spectrum P(k,t) is defined by the expectation value 

(<5(k, t)5*{k', t)) = (27r) 3 5(k - k')P{k, t), (5) 

where the angle brackets denote an average over all realizations. In linear theory, the spatial 
and time dependence of 5 can be factorized, so S(w,t) = 5 (w)D(t)/ D(t ), where t is 
the age of the Universe, 5o( w ) = ^( w ^o)) an d the growth factor (written as a function of 
redshift) is [T3 



5n E(z) too i + z > 

D{z) = —2—J z jm? dz - (6) 



The time t is related to the redshift by 
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t( ) ^ [°° 

Hq J z 



dz' 



l + z')E(z' 



(7) 



Therefore, the power spectrum is P(k,t) = P(k,0)(D/Do) 2 , where D$ = D(to). In an 
Einstein-de Sitter Universe, a/ao = (t/t ) 2 ^ 3 , and D(t)/D = a(t)/a . 

In linear theory, the Fourier components of the velocity field are related to those of the 
density field by 



ia(t) D , . ia(t) D , ~ . 



(8) 



It can also be shown 1113 that 



[l + z] 



Dad 5f2o d 
D~ d~a + ~2Ta[E(z)} 2 D(z) 



(9) 



III. THE OSTRIKER-VISHNIAC EFFECT 

The Ostriker-Vishniac effect is the CMB anisotropy produced by scattering from ionized 
regions or clouds with bulk peculiar velocities. The fractional temperature perturbation in 
the direction 9 induced by bulk motions is 

AT rno . 

p(6) = — (9) = - / n e a T e- T \9 ■ v(w9; w)/c]a{w)dw, (10) 
T Jo 

where n e is the electron density along the line of sight, v(w; w) is the bulk velocity at position 
w at a conformal time 770 — w, <jt is the Thomson cross section, and r is the optical depth 
from us to w. The factor of a(x) arises because the physical time is dt = a(x)dx. [Note that 
9 represents a three-dimensional unit vector along the line of sight, whereas 9 will refer to 
a dimensionless two-dimensional vector in the plane perpendicular to the line of sight — i.e., 
for directions near z, 9 = (0 X , 9 2 , - 9\ - 6%) ~ {9 U 9 2 ,1), whereas 9 = (0i,0 2 ,O).] The 
visibility function, 

g ( w ) = w e (w)a T a(w)e~ T{ - w) = (dr/dw)e- T (11) 

is the probability distribution for scattering from reionized electrons. Here, n^(w) = 
QbPcXe{w)(l + z) 3 /m p is the mean electron density, where f2 6 is the baryon density today in 
units of the critical density (actually, this number should be multiplied by 7/8 to account 
for the neutrons bound in helium nuclei); p c is the critical density; x e (w) is the ionization 
fraction; and m p is the proton mass. The visibility function is normalized so that 

POO 

/ g(w)dw = I - e~ Tr , (12) 
Jo 

where r r is the optical depth to the standard-recombination surface of last scatter, the 
fraction of CMB photons which scattered after standard recombination at z ~ 1100. 
With these developments, the fractional temperature perturbation becomes 
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p(0) = - dwg(w)e-q(w6,w), (13) 



o 



where q(w, w) = [1 + 5(w, w)]v(w, w). According to Eq. ([13]) , only the component of v 
along the line of sight contributes to the anisotropy. In linear theory, we neglect 5^1 
and the total anisotropy is the sum of the contributions of each Fourier component of the 
velocity field. However, in linear theory (and in fact to all orders), v(k) oc k, so only k 
modes along the line of sight can contribute to the anisotropy. The contributions of troughs 
and crests of each Fourier mode cancel approximately when projected along the line of sight. 
Therefore, the contribution of bulk velocities to anisotropies vanishes (or is very small) to 
lowest order in perturbation theory on small angular scales. On larger angular scales, there 
are fewer crests and troughs, so the cancellation is not as complete. 

In other words, angular correlations (on small angular scales) can be due only to k modes 
perpendicular to the line of sight, and these contribute nothing to the anisotropy from bulk 
velocities alone. Furthermore, no anisotropy can be produced by q modes parallel to k, 
and therefore the anisotropy can be due only to the components q_i_(k) perpendicular to k 

lit 

In the Appendix, we make this argument more precise and show that the projection, 
Eq. (|13"1), has Fourier coefficients that satisfy 

(p(«0pV)> = {^?5{R - k')P p (k) } (14) 

with 



1 f^o n 2 (w) 

W = o / P±(k/w, w) dw } (15) 

2 Jo w z 



in a flat Universe [making the replacement Eq. (^) for a closed or open Universe] where 
P±(k,w) is the three-dimensional power spectrum for q_i_(k). We now need to determine 
this three-dimensional power spectrum. 

Again, v(k) oc k, so q_i_(k) can only come from the Fourier transform of the nonlinear 
term <5(w)v(w). The Fourier transform of 5(w, t)v(w, t) [not the full (1 + <5)v] is 

1 r r/ 3 k' r ~ 

q(k, t) = -J — [6(U, t) v(k - k', t) + 5(k - k', t) v(k', t) 

iaDD r d 3 k' ~ w ~ , 1/x /k-k' k'\ 

So(k'Mk-k') [—— + — ). (16) 



2D 2 Q J (2vr) 3 uv ' uv ' V|k-k'| 2 k' 2 
The components of q perpendicular to k are 

q±,i(Kt) = [ $ij - <?j(k,t) 

d3w 5 o( k%(k - k ') ( k - i ( _ _ ___ i . (ir) 





(2tt) 3 uv ' uv ' \ 1 k 2 J \k' 2 |k-k 
For a Gaussian density field 



6 



(5o(ki - k%(k')~5*(k 2 - k>*(k")> 

= (27r) 6 P(|ki - k'|,t)P(A;',t)[(5(k 1 - k 2 )5(k' - k") + 5(ki - k 2 )5(ki - k' - k")]. (18) 
We combine these expressions to find (after some extensive although straightforward algebra) 
(q ± (k, t) • ql(k', t)) = (2tt) 3 5(k - k')P ± (k, t), (19) 

with 

PxM = ^(§)V); (20) 

the Vishniac power spectrum is || 

S(k) = kfdy £ dxP(ky)P(kyfl + y 2 - 2xy) ( \~ fj^f ■ (21) 

(The differences with Eq. (2.13) in are the result of differing Fourier conventions.) 

Assembling these results, in a flat Universe, the Ostriker- Vishniac effect produces a power 
spectrum, 

PM = JL jT ^ [a(„)] 2 ^ S(k/w) dw, (22) 

with the Vishniac power spectrum S(k) given by Eq. fl2iP . In an open or closed Universe, 
one makes the replacement Eq. (|52| ) in the argument of S(k/w) and in the two factors of w 
in the denominator of the integrand. the numbers for the open models in the table and the 
open-model curves need to be fixed. 



A. Gaussian Visibility Function 



If the visibility function is approximated as a Gaussian in conformal time, 



exp 



1 [w 

2 (Aw) 2 



(23) 



where w r is the comoving distance to the re-scattering surface, Aw is the width of the re- 
scattering surface, and r r is the fraction of CMB photons re-scattered, then the CMB power 
spectrum is 



P p {k) 



1-e 



327r 2 v /7r(Au>)u> 



a(w r )D(w r )D(w r ) 



n 2 



S(n/w r 



(24) 



In practice, the visibility functions for the ionization histories usually considered are very 
poorly approximated by a Gaussian. Still, the Gaussian approximation provides a simple 
analytic description of how the anisotropy spectrum depends on the redshift of re-scattering, 
optical depth, and S(k). 
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B. Relation to Multipole Moments 



Given the power spectrum P p (k) and our Fourier conventions, it is straightforward to 
write the temperature autocorrelation function in terms of the power spectrum: 

C(a) = (p(0 + a)p0j) = — I ndn J q (ko) P p {k). (25) 

27T J 

The correlation function can also be written in terms of the more commonly used CMB 
multipole moments Ce as 

2/ 4- 1 

C(a)= y £-^—P i (cosa)C e , (26) 

where P^cosa) are Legendre polynomials. Since P^(cos#) — > Jq(£9) as t — > oo, it follows 
that 

C t = P p (k = £). (27) 

This identification allows us to illustrate our results in the multipole moments that have 
become familiar and to compare our results with those of previous authors. 



C. Comparison with Previous Results 

First, we compare with the results of Ref. |7| (and therefore Ref. ||). To do so, we 
consider only the Einstein-de Sitter Universe as they do. In an Einstein-de Sitter Universe, 
Eq. ( p2"D becomes (after changing the integration variable to rj = r] — w), 

r,(«) = h f ^ W = V^ (" V s ( —) **- ^ 

4vr 2 7o (vo~v) 2 Vo \VoJ \Vo ~ V J 
We then note that they approximate the power spectrum by evaluating S(k/w) at w = r] 
and taking it outside the integral in Eq. fl2"8"D- Now we compare our result with Eq. (6.15) 
in Ref. and their J in their Eq. (6.18). If we note that their K(k)P 2 (k) is equal to our 
S(k)/k and assume that i] — rj ~ i]q in the integrand, then our results reproduce theirs 
(although it appears that their result is twice ours; we have not been able to locate the 
source of the discrepancy). 

Now we want to compare our results with JTO]. Although their derivation, unlike ours, 
relies on a solution to the Boltzmann equation, our results should agree. Specifically, our 
Eq. Q2"2"D should agree with their Eq. (24). To check that this is so, we change our integration 
variable in Eq. (|22|) to k = £/w and find 



C t = I dkW e (k)S(k) } (29) 

where 



'A) 



/ \ 2 

w ^ = T^ [9iw = i/k)]2 lw) ■ (30) 

If we put their Eq. (24) into this form [to do this, keep in mind that their overdot denotes 
differentiation with respect to r\ instead of t, so our aD is equal to their D — also, they 
normalize D such that D(to) = 1], we find agreement with their results. 
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D. Variance of Temperature Distribution 



The variance of the temperature distribution induced by the Ostriker-Vishniac effect is 



((AT/T) 2 ) = -L J £d£C e . 
In the Gaussian approximation for g(w) and in a flat Universe, 



((AT/T) 2 ) 



a(w r )D(w r )D(w r 



327r 3 y7r(Au;) 

and for an Einstein-de Sitter Universe this simplifies to 



((AT/T) 2 ) 



t.: 



4 fv 



(31) 



k dk S(k), 



(32) 



)0 
/ kdkS{k). 



(33) 



IV. THE SPATIAL POWER SPECTRUM 



For the power spectrum, we use 



P(k) 



~8 



S 2 H (k/2) n T 2 (k p Mpc/hT), 



(34) 



where T(q) is the usual CDM transfer function, k p = ka^ = kH /2 is the physical wave 
number with our conventions, and T ~ QqIi is given more accurately in terms of Qoh and 
the baryon fraction Qj, by Eqs. (D-28) and (E-12) in Ref. The factor of 8 in the 

denominator in Eq. (0) arises because we are using a Ho = 2. For the transfer function, 
we use rirOT 



T(q) 



ln(l + 2.34g)/(2.34g) 



[1 + 3.89g + (16.L?) 2 + (5.46g) 3 + (6.71g) 4 ]V4 • 



(35) 



For 5h, we use the fits to the COBE anisotropy given in ||T^]. For the tilted model, we use 
the COBE normalization from Ref. Jl7j obtained including tensor modes with the amplitude 
predicted by power-law inflation. (Note that the tensor contribution in other inflationary 
models may be different.) 

Alternatively, the power spectrum may be normalized at small distance scales through 
the cluster abundance which fixes a§, the variance in the mass enclosed in spheres of radius 
8 hr 1 Mpc, to cr 8 ~ (0.6 ± 0.1)f2o °' 53 H - In terms of the power spectrum, 



2tt 2 



k 2 dkP(k) 



3ji(k p R) 



(36) 



where R = 8h 1 Mpc, and ji(x) is a spherical Bessel function. Since we are using a Q ^ 1, 
k p (rather than k) enters into the argument of the spherical Bessel function. 
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V. NUMERICAL RESULTS 



In this Section we present numerical results for the anisotropy spectrum and rms temper- 
ature fluctuation for a variety of CDM models and ionization histories. We consider a family 
of ionization histories parameterized by a constant ionization fraction x e after a redshift z r . 
In a cosmological-constant Universe, this induces an optical depth, 



r r = {yt b p c a T a Q /m, 



'l + z) 2 dz 



0.069 Q b hx e r ^ + J^ dZ . (37) 



In a flat Universe, this becomes 

r r = 0.046 Qbhx e 



o E(z) 

(1 + ^) S 
o E(z) 



fi (l + Zr) 3 + (1 " n ) ~ 1 



(38) 



The prefactor 0.069 (or 0.046 in the flat-Universe expression) is obtained if all the baryons 
are in the form of protons, an assumption often made when calculating ionization histories. 
If one takes into account the fact that one quarter of the baryonic mass is helium, then the 
prefactor should be multiplied by 7/8. 

Table | lists the predicted rms temperature anisotropy ((AT) 2 ) 1 ^ 2 due to the Ostriker- 
Vishniac effect for a variety of CDM models, which are parameterized by Qq, h, n, and 
Qbh 2 , and ionization histories, which are parameterized by x e and z r . The Table also lists 
r r , ag, and Ugfig 53 . We present a variety of models and ionization histories to illustrate the 
effects of variation of various parameters on the predictions. It should be kept in mind that 
not all of the CO-BE-normalized models satisfy the aforementioned observational constraint 
to a 8 from cluster abundances [IB|. We therefore list in the Table the predicted RMS tem- 



perature anisotropy expected if the COBE normalization is disregarded and the amplitudes 
of the power spectrum is fixed instead to the cluster abundance. Similarly, the h = 0.35 
model is meant to be illustrative; this value of the Hubble constant is in disagreement with 
most measurements. 

Figure [I] shows the predicted anisotropy for the canonical standard-CDM model (fl = 1, 
= 0, h = 0.5, n = 1, Q^h 2 = 0.0125), for a variety of ionization histories, with the 
linear power spectrum, also shown, normalized to COBE. We know from quasar absorption 
spectra that the Universe has been significantly reionized at least since a redshift of z r = 5. 
Therefore, we include predictions for this minimal level of anisotropy. 

The effect of reionization on the primary anisotropics can be quantified primarily by the 
optical depth r r [Q; any two ionization histories that have the same r r produce roughly the 
same primary anisotropies. However, as Fig. [I] and Table | show, the secondary anisotropies 
are more dependent on the ionization fraction than on the optical depth. 

The growth of density perturbations leads to a growth of the peculiar velocities that 
induce the secondary anisotropy. Therefore, if reionization takes place later, the anisotropy 
is significantly larger, even for the same optical depth. Furthermore, if reionization takes 
place later, the peak of the anisotropy spectrum, which is fixed in comoving distance scale 
by the peak of the power spectrum, subtends a larger angle. In open models, the peak of the 
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TABLE I. The rms temperature perturbation ((AT/T) 2 ) in //K due to the Ostriker-Vishniac 
effect and value of £ pea k at which t{l + l)Cg peaks for a variety of CDM models and ionization 
histories for models normalized to both COBE and cluster abundances. 
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FIG. 1. Multipole moments for the Ostriker-Vishniac effect for the COfiE-normalized canon- 
ical standard-CDM model (ft = 1, h = 0.5, n = 1, ft b /i 2 = 0.0125), for a variety of ionization 
histories, as listed. We also show predictions for several open high-baryon-density models with the 
same x e and r r , normalized to the cluster abundance, with dashed curves. The dotted curves show 
the primary anisotropy for this model for r r = 0.0, 0.1, 0.5, 1, and 2, from top to bottom. 
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power spectrum moves to larger physical scales, while the effects of geometery shift these 
physical scales back to smaller angular scales; the result is that the spectrum peaks at only 
slightly lower I than in a flat universe. 

The weakness of the dependence of the anisotropy amplitude on r r indicates that a 
good fraction of the anisotropy is produced at later times when the amplitude of the power 
spectrum has grown. The convergence of the predictions at small I for ionization histories 
with the same x e indicates that, in particular, the anisotropy at larger scales comes primarily 
from lower redshifts [pi , and therefore that assuming that re-scattering occurs primarily 
at the epoch of reionization |],[JJ does not provide a good approximation, especially at 
larger scales. Moreover, in open models the evolution of the growth factor differs from a 
flat Universe by larger amounts at later times; therefore the slope of Ce at larger scales is 
greater than in a flat universe. 

The growth of density perturbations induces an initially Gaussian distribution of pertur- 
bations to become non-Gaussian. Therefore, the distribution of anisotropics produced by 
the Ostriker-Vishniac effect will be non-Gaussian even for Gaussian initial conditions, and 
the higher n-pt temperature correlation functions may be calculated in perturbation theory 

The anisotropy is due primarily to density perturbations on distance scales closer to those 
probed by galaxy surveys than to those probed by COBE. Furthermore, the anisotropy is 
due to a large extent to scattering at later times. Therefore, when the power spectra are 
normalized to the cluster abundance, the spread in predicted values for the rms anisotropy 
is smaller than it is if the models are normalized to COBE. Since the anisotropy is a second- 
order effect, the rms anisotropy is proportional to erf. In a flat Universe, models normalized 
to cluster abundances produce rms temperature anisotropies listed in Table [I] of 0.8-3 /iK 
for a constant ionization fraction of unity, whereas an ionization fraction of 0.2 yields rms 
anisotropies of 0.3-0.8 /zK. 

Our numerical results can also be easily scaled for a different value of the ionization frac- 
tion x e if the ionization epoch is held fixed and if r r <C 1. The rms temperature anisotropy 
should be roughly proportional to x e if z r is held fixed. Also, if we neglect the effect of Q b and 
h on the power spectrum (in particular, Q b should have only a relatively weak effect on the 
power spectrum as long as Q b <S ^o)> then Eq. (^) shows that the rms anisotropy should be 
proportional to the combination Q b hx e = (Q b h 2 )x e /h, itself proportional to the column den- 
sity of electrons to the last scattering surface. This also suggests that the Ostriker-Vishniac 
anisotropy in models with a high baryon density should be larger. To illustrate (and to 
illustrate the effects of geometry), we also include in Table | and in Fig. p] predictions for 
an open model with fl b ~ O /3, as suggested possibly by a peak in the observed power 
spectrum at 100 h~ l Mpc |[20|| . High-baryon-density models will be explored more carefully 



in Ref. 21 



VI. APPROXIMATE VISHNIAC POWER SPECTRA FOR CDM MODELS 

Here we provide simple approximate analytic fits to the Vishniac power spectrum S(k) for 
standard-CDM power spectra. These can be used to quickly estimate the Ostriker-Vishniac 
anisotropy for CDM models with an arbitrary ionization history. 
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With the standard-CDM power spectrum given in the form of Eq. fl3~4|), S(k) may be 
written, 

S(k) = ^ (k/2) 2n kI n (k p Mpc/hT), (39) 
16 



where 



IrM) = / q °° dy £ dx(yy/l +y>- 2xy) n T 2 {yq)T 2 {q^ I + y* - 2xy) ^? (4Q) 

For the standard-CDM transfer function, Eq. fl35|), the functions I n {q) may be approximated, 
for 0.5 ^5 n ^ 1, by an asymmetric Gaussian: 



logiofe] -log 10 [g (n)] n " 



(41) 



with log 10 [/i(n)] = -4.716 - 0.636(n - 1) + 0.832(n - l) 2 , log 10 [? (n)] = 0.477 + 1.36(ra - 
1) + 0.975(n - l) 2 , and a = 0.6 + 0.33(n - 1) + 0.34(n - l) 2 for q < q , and a = 0.87 + 
1.105(n - 1) + 1.19(n - l) 2 for q > q . 

Fig. |2| shows the multipole moments obtained from numerical calculation of S(k) and 
the analytic approximation given here for a variety of models. Although this analytic ap- 
proximation does not accurately reproduce the asymptotic large- and small-A; behavior of 
S(k), it accurately approximates the peak of k 2 S(k). As indicated by Eq. (p2|), this should 
provide reliable estimates of the mean-square anisotropy. 



VII. DISCUSSION 

Although precise limits depend on imprecisely determined cosmological parameters, the 
existence of significant CMB anisotropy at the degree scale and current limits to arcminute- 
scale anisotropy and distortions to the CMB blackbody spectrum seem to indicate that 
T r % 1. However, some reionization is generically expected in CDM models. Again, the 
detailed ionization history will depend on the model and on details of star formation that 
are still difficult to quantify reliably. Still, optical depths r r = (9(0.1) are observationally 
consistent and expected in CDM models. 

Such a level of reionization will lead to secondary anisotropies at arcminute scales like 
those shown in Section V. Furthermore, the bulk of the reionization in CDM models generally 
occurs later, so the secondary anisotropy will be larger for a given r r than if reionization 
occurred earlier. The predicted level of anisotropy in virtually all of the models we considered 



here is smaller than the current limit of about 50 fiK at arcminute angular separations |22 
However, there has been much recent progress in the design of interferometer experiments 
dedicated to CMB measurements, and there should be much more in forthcoming years. To 
estimate the signal that might be detectable, we suppose that a detector with a sensitivity s 
(in fiK A/sec) maps a fraction / sky of the sky with an angular resolution 9f w h m < 1' sufficient 
to detect these anisotropies for a time t mon months. Such a detector should be sensitive to 
an rms temperature anisotropy as small as 
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FIG. 2. Semi-analytic fits to S(k), the Vishniac power spectrum. Top: full calculation (solid) 
and the analytic approximation (dashed) for the standard-CDM S(k). Bottom: £(£ + l)C^/(27r) 
calculated with the full S(k) (solid) and the approximation (dashed) for standard CDM with an 
ionization history given by z r = 19, x e = 1.0, r r = 0.1. 
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Although the interferometry experiments that will likely reach these angular scales are not 
really parameterized by a sensitivity s, the effective value of s for these experiments is already 
better than 1000 fiK y/ sec. Furthermore, with a full power-spectrum analysis of the data, the 
sensitivity could be improved perhaps by an order of magnitude over this simple estimate. 
Comparing Eq. (^) for a reasonable t mon with the predictions in Section V, it seems quite 
plausible — even though the experiments may be challenging — that the Ostriker-Vishniac 
anisotropy expected in many leading CDM models will be detectable. The results of Ref. 
|TT| imply that the anisotropy generated by some of the more optimistic models considered 



here may already be detectable by upcoming experiments such as CBI, especially if binned 
over a wide range of £. 

Observation of the Ostriker-Vishniac effect will be important for understanding the ion- 
ization history of the Universe and for determining cosmological parameters with the CMB. 
The linear-theory anisotropies depend primarily on the optical depth r r and only secondarily 
on the details of the ionization history (e.g., the values of z r and x e ). Therefore, degree-scale 
anisotropies may determine r n but they will not strongly constrain the epoch of reionization. 
When combined with sensitive measurements of the degree-scale anisotropy, measurement 
of the Ostriker-Vishniac effect will therefore provide the epoch of reionization. 

The shape of the CMB power spectrum for these anisotropies is relatively featureless. 
Therefore, even detection of the anisotropy, without a precise determination of the shape of 
the power spectrum, will provide the most essential information on the ionization history. 
Of course, more precise constraints to the ionization history can be obtained if the power 
spectrum can be mapped and, some shape information along with multi-frequency obser- 
vations will be necessary to distinguish the Ostriker-Vishniac anisotropy from foreground 
contamination such as that induced by extragalactic point sources. 

Determination of the geometry of the Universe relies primarily on the first acoustic peak 
in the primary CMB spectrum fJTJ,^3|| . However, precise and reliable determination of other 
cosmological parameters, such as the baryon density, Hubble constant, cosmological con- 
stant, and spectral index, depends on the structure of the smaller-scale acoustic peaks. The 
details of the ionization history (e.g., the epoch of reionization, which we have explored here) 
will, in fact, have some effect on these higher peaks. This will in turn affect the reconstruc- 
tion of some cosmological parameters, although this uncertainty has not yet been taken into 
account in the parameter-estimation analyses of Refs. [OL Therefore, the information on 
the epoch of reionization provided by mapping the secondary anisotropies will be important 
for precise determination of these other cosmological parameters with MAP and Planck. 

It seems unlikely that the current FIRAS constraints to Compton-?/ distortions of the 
CMB blackbody spectrum will be improved significantly. Therefore, the Ostriker-Vishniac 
effect seems to be the most promising probe of the ionization history. In addition to its sig- 
nificance for cosmological-parameter determination, the Ostriker-Vishniac effect can provide 
a window to the epoch of the earliest collapsed objects that is inaccessible with any other 
observations. 
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APPENDIX: PROJECTED POWER SPECTRUM FOR A DIVERGENCE-FREE 

VECTOR FIELD 

If an isotropic random field q(w) is decomposed into Fourier components, 

q(k) = J d 3 w exp(— zk • w) q(w), (43) 

then the Fourier components may be written as the sum of divergence-free and curl-free 
parts parts: q(k) = q_i_(k) + qy(k) with k • q_i_(k) = and k x q|(k) = 0; in particular, 
q||(k) = kg||(k). 

Suppose we observe a projection 

-> rvo 

p(9) = — dw g(w) 9 ■ q(w6,w;w), (44) 

J 

of 9 ■ q(w9, w), the component of q(w) along the line of sight, with a visibility function g(w), 
and we define the Fourier components, 

p(R) = J d 2 9p{9)e- lfiS . (45) 

The contribution of the parallel components q(k) to p(9) will be small because crests and 
troughs of each Fourier component projected along the line of sight will tend to cancel. 
The purpose of this Appendix is to show that if the perpendicular components satisfy 
(q_i_(k) ■ qj_(k')) = (27r) 3 5(k — \t')P±(k), then the Fourier components of the projection 
satisfy 

(P(k)P*(k')) = (2tc) 2 5(k - k')P p (k). (46) 

with 



1 rvo a 2 (w) 

P p {k) = - / dw g -\^P^K/w,w). (47) 
— J xv 



To do so, we follow the steps in the Appendix of Ref. p4| that lead to the Fourier-space 



analog of Limber's equation. Rather than reconstruct the entire argument, we refer the 
reader to Kaiser's paper [^3] for more details. 



The components of an isotropic divergence-free field satisfy 

(q ±>i (k)T ±jj (k')) = hSq - k^-XcUk) ■ ql(k')). (48) 
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If the contribution of a narrow shell of width Aw; centered at wo is 



pwq+Aw/2 

Ap(6) = g(w ) / dw9 ■ q±(w 9 1 ,w 9 2 ,w;w), (49) 

J wn—Aw/2 



rWQ+Aw/2 
I wo— Aw/2 

then with Eq. (|48|) , we can show that its Fourier coefficients Ap(H) satisfy 

(Ap(k) Ap(H')) = (2tt) 2 5(k - k') AP p {k). (50) 

with 



APp{K)= ii^ [i-(«. t)i[ fc (»^./ 2 )]'i. ± y5 + «), (si, 

and k = (ki/w , ^2/^0, k 3 ). As argued by Kaiser, only modes with k 3 <C n/w (those 
very nearly perpendicular to the line of sight) will contribute appreciably to the integral, so 
(9 ■ k) 2 <C 1. We then add the contributions Ap(/c) along the line of sight to obtain Eq. (pETp . 
In an open or closed Universe one replaces 



S(a H wJ\l - n - n A \) 
w — > / =— (52) 

O'qHqJ 1 — VLq — 



in the first argument of P(k/w,w) in Eq. (f47|), and in the two factors of w that appear in 
the denominator, where S(x) = sinhx in an open Universe and S(x) = sinx in a closed 
Universe 0.0 

Note that this derivation also explains why the contribution from qn is suppressed. This 
term satisfies 



(qyik^ik')) = kkW^')) = (2vr) 3 ^(k - k%-^|] (k). (53) 

This would therefore have a contribution to the projected power spectrum of ^-q proportional 
to(£-k) 2 <l. 

The approximation used here is valid for small angular separations, k ^> 1. Since the 
Ostriker-Vishniac effect produces anisotropies at multipole moments i > 1000, this flat-sky 
approximation should lead to errors no greater than 0(0.1%). At large angular scales, we can 
also expect the anisotropy from the parallel term to be significant [pTojl . Also, the visibility 



function g(w) should be smoothly varying over wavelengths at which there is significant 
power. Again, the visibility function for re-scattered CMB photons will be very broad, so 
this approximation should be fine for our calculation. 



1 We thank Wayne Hu for pointing out an error in this scaling in a previous draft. 
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